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Discrete time random walks, in which a step of random sign but constant length Sx is performed 
after each time interval St, are widely used models for stochastic processes. In the case of a correlated 
random walk, the next step has the same sign as the previous one with a probability q^\. We extend 
this model to an inhomogeneous ensemble of random walkers with a given distribution of persistence 
probabilites p(q) and show that remarkable statistical properties can result from this inhomogenity: 
Depending on the distribution p(q), we find that the probability density p(Ax, At) for a displacement 
Ax after lagtime At can have a leptocurtic shape and that mean squared displacements can increase 
approximately like a fractional powerlaw with At. For the special case of persistence parameters 
distributed equally in the full range q£ [0, 1], the mean squared displacement is derived analytically. 
The model is further extended by allowing different step lengths Sxj for each member j of the 
ensemble. We show that two ensembles [St, {{qj, Sxj)}] and [St 1 , {{q'j, Sx'j)}] defined at different time 
intervals St ^ St' can have the same statistical properties at long lagtimes At, if their parameters 
are related by a certain scaling transformation. Finally, we argue that similar statistical properties 
are expected for homogeneous ensembles, in which the parameters (qj(t),Sxj(t)) of each individual 
walker fluctuate temporarily, provided the parameters can be considered constant for time periods 
T 3> At longer than the considered lagtime At. Similar models are applicable to many complex 
systems in which the individual agents undergo distinct - yet aysnchronous - behavioural phases, so 
that the statistics of the ensemble as a whole can still be considered as stationary. 
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I. INTRODUCTION 

Many natural phenomena can be described as continu- 
ous stochastic processes. The prototypical example for a 
stochastic phenomenon is Brownian diffusion of a parti- 
cle within a liquid, usually modeled as a Wiener process. 
Other systems, such as a particle diffusing within a ho- 
mogeneous gravitational field, involve a combination of 
stochastic and deterministic forces. These cases can be 
modelled by adding a drift term to the stochastic dif- 
ferential equation of the Wiener process. Furthermore, 
many systems show a tendency to return to a preferred 
equilibrium state, such as in the case of particle diffusion 
within a harmonic potential well, which is modeled by 
the Ornstein— Uhlenbeck process process. 

Stochastic processes such as the Wiener or 
Ornstein— Uhlenbeck process deal with a continuous 
state variable (e.g., the particle's coordinate x(t)) 
as a function of continuous time. However, when a 
continuous particle trajectory is measured, the re- 
sulting data is typically a time series Xt, where the 
index t G {0,1,2,...} is counting discrete time points 
separated by a fixed sample interval 5t. It is most 
natural to model such time series as a discrete stochastic 
processes, also called a Discrete Time Random Walk 
(DTRW) model. While, in general, a DTRW may have a 
continuous state variable x G 5R, it is often advantageous 
to further simplify the mathematics and to reduce also 
the state space of the system to a discrete set, such as 
x t G {. . . , —2, —1, 0, 1,2,.. .}. The simplest example is 
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an uncorrelated random walk where the state variable 
x t changes after each time interval St by plus or minus 
a fixed step length 61. This model can be shown to 
converge towards the Wiener process in the limit of 
small St and small 51. 

Possible applications of DTRW models go far beyond 
traditional physics and include, for example, many bio- 
logical and economic systems. In particular, the spatio- 
temporal motion patterns of various biological agents 
(animals, cells) could be successfully described by DTRW 
models, but only after correlations between successive 
steps have been included. In the simplest version of a dis- 
crete time correlated random walk (DTCRW), the next 
step has the same sign as the previous one with a prob- 
ability q 7^ |. For values q < \, the random walk is 
called anti-persistent, correponding to a back-and-forth 
motion. For q = 4 it behaves like an uncorrelated ran- 
dom walk. For q > | , the random walk is called persis- 
tent, since it consists of long chains of steps in the same 
direction. This model is equivalent to a Markov chain 
and has first been investigated in 1920 [T]. It is possible 
to derive analytically the resulting displacement distribu- 
tion P(Ax, At), which is defined as the probablity that 
the 'particle' has moved a distance of Ax — kSx during 
a time interval At = n8t, as well as the mean squared 
displacement Ax 2 (At) [2J. 

Besides Sx and St, which set the spatial and tempo- 
ral scale of the system, the only true parameter of this 
standard correlated random walk model is the persistence 
parameter q G [0,1]. It directly determines the correla- 
tion time t = 2(i- q) °f t ne velocity autocorrelation func- 
tion C vv (At) = <7 2 e~ A */ T . This correlation time r be- 
comes also apparent in the mean squared displacement 
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Ax 2 (At) as the characteristic time separating the bal- 
listic regime (small lagtimes) from the diffusive regime 
(large lagtimes). 



II. MODEL AND RESULTS 

Standard Correlated Random Walk: CRW 
model 



(a) 




FIG. 1. Example of a temporally inhomogeneous random 
walk (b), consisting of three subsequent correlated random 
walks with different degrees of persistence (a). 



Not all correlated random walks, however, can be char- 
acterized by a single correlation time r. For example, 
it is known that individual cells that are migrating on 
planar substrates can switch between different migration 
modes, ranging from anti-persistent to highly persistent 
motion [3], [4]. The mean squared displacement of such 
cells often approximates a fractional powerlaw over sev- 
eral decades of lagtime [5], pointing to a scale-free ran- 
dom process without any unique correlation time. In or- 
der to describe such temporally inhomogeneous systems, 
it is natural to use so-called superstatistical models [5] , in 
which the statistical parameters (correlation times, dif- 
fusion constants, etc.) are themselves subject to random 
temporal fluctuations (compare Fig{T]). 

In this paper, we start from the standard correlated 
random walk (CRW) model with persistence param- 
eter q and its known quantities P(Ax, At \ q) and 
Ax 2 (At | q). We then generalize the model by allow- 
ing q to vary throughout the ensemble, according to a 
fixed distribution P(q). The resulting displacement dis- 
tribution (P(Ax, At)) q and mean squared displacement 

(Ax 2 (At)) of the ensemble can simply be obtained 



by computing the average of these q-dependent quanti- 
ties, weighted with P(q). The average will be performed 
numerically, revealing almost exponentially shaped dis- 
placement distributions and a mean squared displace- 
ment resembling a fractional powerlaw over several orders 
of magnitude in the lagtime. Furthermore, we argue that 
similar statistical properties are expected for homoge- 
neous ensembles, in which the parameters (qj(t),Sxj(t)) 
of each individual walker fluctuate temporarily, provided 
the parameters can be considered constant for time peri- 
ods T 3> At longer than the considered lagtime At. 



Since the statistical properties for this model have been 
derived analytically in Ref. [2] , we just repeat the main re- 
sults here. The displacement distribution for a DTCRW 
with persistence parameter q is given by 



P(k,n\q) = P(Ax = kSx, At = nSt | q) 

{n-\k\)/2 
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where n and k must either both be even or both be odd. 
It has additionally been assumed that initially (at time 
step t = 0) the probablity for the particle to go left or 
right are equal. Note that P(k, n) — for \k\ > n. 
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FIG. 2. Displacement distributions in the correlated random 
walk model. Case (a) correponds to an anti-persistent walk 
(q = 0.1), case (b) to the uncorrelated walk (q — 0.5) and 
(c) to a persistent walk (q — 0.9). In all cases, the distribu- 
tions are shown for the three lagtimes n = 8 (circles), n = 64 
(squares), and n = 512 (stars). The distribution is only de- 
fined at the discrete points, lines are guides for the eye. Note 
that P(k,n) = for \k\ > n. 

Fig |2] shows the evolution of the displacement distribu- 
tion with lagtime for an anti-persistent, the uncorrelated 
and a persistent case. For lagtimes much larger then the 
correlation time, the distributions approach Gaussians. 
The mean squared displacement is given by 



Ax 2 (n | q) 
n q 



8x 



Ax 2 (At = nSt | q) = 
(2g-l)[l-(2g-l) 



(2) 



2ng(l - q) 

Fig{3] shows the mean squared displacement (MSD) as 
a function of lagtime for the same three degrees of persis- 
tence as in Fig[2j In the uncorrelated case (green line), 
the MSD is linear, corresponding to diffusive random mo- 
tion. An anti-persistent walk (red line) with the same 
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lagtime n 

FIG. 3. Mean squared displacement in the correlated random 
walk model. The red line correponds to an anti-persistent 
walk (q = 0.1), the green line to the uncorrelated walk (q = 
0.5) and the blue line to a persistent walk (q = 0.9), just as 
in Fig|2] The thin black lines show a linear and a quadratic 
dependence. 



step length has a smaller MSD, compared to the diffu- 
sive case. For a persistent random walk (blue line), the 
MSD starts ballistically for lagtimes smaller than the cor- 
relation time and then continues to grow diffusively. 



culated analytically. For example, in the case of an equal 
distribution P(q) — 1 in the full range q£ [0, 1], the av- 
eraged mean squared distribution is given by 



Ax 2 (n | q) 



e+(n + l) 



q.P{q) = l 

r , n + 3 — e 



(5) 



Here, e = for even times n and e=l for odd times n. 
The digamma function ^(x) — 4- log (T(x)) is defined as 
the logarithmic derivative of the gamma function. For a 
derivation of Eq. ^ see Appendix. 
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FIG. 4. Displacement distributions for an inhomogeneous 
ensemble of correlated random walkers. Case (a) corresponds 
to equally distributed degrees of persistence with P(q) = 1, 
case (b) to a distribution P(q) oc j 4 . In all cases, the dis- 
tributions are shown for the three lagtimes n = 8 (circles), 
n — 64 (squares), and n — 512 (stars). The distribution is 
only denned at the discrete points, lines are guides for the 
eye. Note that P(k,n) = for |fe| > n. 



B. Inhomogeneous Ensemble 

We now assume an inhomogeneous ensemble of cor- 
related random walkers j with different values qj, dis- 
tributed according to a given distribution p(q). We can 
then compute the statistical quantities of the ensemble by 
averaging over the corresponding quantities of the stan- 
dard correlated walk, so that 

(P(k,n)) = f dqP(q)P(k,n\q) (3) 
Jg=0 

and 



We have also investigated distribution of the form 
P(q) oc q 1 for the persistence parameter q G [0, 1] and 
computed the averages numerically for the cases 7 = 
and 7 = 4. The resulting distributions (Figj4]) are 
strongly leptocurtic and resemble exponential functions, 
rather than Gaussians, which is usually interpreted as a 
signal of anomaleous random motion. 

The corresponding mean squared displacements 
(Fig{5]) also show clear signatures of anomaleous behav- 
ior. In particular, for the case of equally distributed de- 
grees of persistence (blue line) , the MSD can be well ap- 
proximated by a powerlaw with a fractional exponent of 
about 1.3 (middle dashed line). In order to confirm our 
numerical averaging procedure, we have also compared 
this MSD curve with the analytical result from Eq.([5} 
and found excellent agreement (black dots). 



(Ax 2 (n I q)) = [ dq P(q) Ax 2 {n | q). (4) 

^ ' 1 Jg=0 

There are specific choices for the distribution P(q) of 
peristence parameters for which the integrals can be cal- 



C. Extension of the model and scaling poperties 

A natural extension of the above model is an inho- 
mogeneous ensemble where each walker has not only its 
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FIG. 5. Mean squared displacement for an inhomogeneous 
ensemble of correlated random walkers, corresponding to the 
distributions P(q) oc q 1 of Fig|4] The blue line despicts the 
case 7 = 0, the green line the case 7 = 4. The dots show the 
analytical result according to formula Eq.|5|. The thin black 
lines show the powerlaws n 1 , n 1 ' 3 and n 2 . 

own persistence parameter qj , but also an individual step 
length Sxj . This extension is actually required for trans- 
forming a given ensemble of random walkers from sam- 
pling interval St to another sampling interval S' t , so that 
the statistical properties (such as P(Ax, At \ q) and 
Ax 2 (At | q)) remain invariant: It has been shown in [7] 
that a single DTCRW with parameters [St, q, dx] has the 
same statistical properties as a DTCRW with rescaled 
parameters [St' = St/s, q' ,Sx'], if the following scaling re- 
lations are used: 

g' = g , (g,a) = i(l±|2«-l| 1 /*),and (6) 



Sx = Sx ■ g(q, s) with g(q, s) = J -jyz %r-- (7) 

V IK 1 ~ 1) s 

These scaling relations have to be applied to each in- 
dividual member of the ensemble. When transforming 
the ensemble to a smaller sampling time St' < St, all 
step sizes decrease, the anti-persistent walkers become 
even more anti-persistent and the persistent walkers even 
more persistent (compare Figj6f. We have tested the 
correctness of the scaling transformation by comparing 
the displacement distributions of an ensemble sampled 
at St = 1 and the corresponding rescaled ensemble sam- 
pled at St = 1/8 (compare Fig{7]). 
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FIG. 6. Scaling transformation of an inhomogeneous en- 
semble of correlated random walkers to different sampling in- 
tervals 8t, using the scaling transformations Eq[5] and Eq[7| 
Each individual member j of the ensemble is represented as a 
point (qj,Sxj). Solid lines are guides for the eye. For St = 1, 
the persistence parameters qj are equally distributed in [0,1] 
and all step lengths Sxj are set to 1. 
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FIG. 7. Displacement distributions of an inhomogeneous 
ensemble of correlated random walkers for three different lag- 
times. Solid lines correspond to the ensemble sampled at time 
increments of St — 1. The dots show that the distributions 
remain invariant when the ensemble is rescaled to a smaller 
sampling time St = 1/8 according to the scaling transforma- 
tions Eq|6] and Eqj7] 
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D. Ensembles of temporally inhomogeneous 
walkers 



Above, we have considered an inhomogeneous ensem- 
ble of random walkers with temporally constant prop- 
erties. However, similar statistical properties would be 
expected in the case of a homogeneous ensemble of ran- 
dom walkers with temporally fluctuating properties. For 
a concrete example, consider a Petri dish with a set of 
migrating cells, all of the same type and prepared in the 
same way. Assume that each specific cell j of this ensem- 
ble switches its persistence parameter q in regular inter- 
vals T = NSt according to the distribution p(q) and pro- 

duces a time-series x t corresponding to the x-coordinate 

(i) 

of its trajectory. Such an individual time-series x t does 
not represent a stationary random process, since the sta- 
tistical properties (such as the displacement probability 

P(xt2m~ x t )) are not time-invariant but depend on the 
momentary q- value of cell j. However, if the individual 
cells switch their q- values in a completely independent, 
asynchronous way, we can compute the ensemble-average 
(/p(s^ m — x^)^ ) and thus obtain a stationary quan- 
tity. The q- average used above is just a convenient way 
to compute this ensemble-average. 



III. DISCUSSION AND OUTLOOK 

An important point to be addressed is the nature of 
the approximate 'powerlaw' found in specific cases such 
as shown in Figj5] It is clear that our random process 
is not really long-time-correlated and scale-invariant (not 
even after performing the ensemble-average) since all cor- 
relations are lost for lagtimes larger than the duration T 
of each constant-g-phasc. Nevertheless, this duration T 
can last over several decades of lagtime and may exceed 
the maximum measured time interval in an experiment. 

In this paper, we have considered just a very simple 
case of a superstatistical random walk. The basic concept 
of starting with a homogeneous random walk model, con- 
trolled by a set of parameters qi , q 2 , . . . , and then defining 
a second layer of random 'super-process' for the dynamics 
of these control parameters, leads to a multitude of pos- 
sible model variants. These variants may be particularly 
useful for biological systems with extremely inhomoge- 
neous ensembles. An interesting problem will be how to 
extract from a measured superstatistical time-series the 
properties of the super-process from those of the under- 
lying random walk. 



IV. APPENDIX 

We consider a DTCRW model with 8x = 1, St = 1 and 
persistence parameter q. The velocity autocorrelation 



function is given by 



C vv {n I q) = = (2q I)™. (8) 



For the case of equally distributed degrees of persistence, 
p(q) = 1, it is straight forward to calculate the average 
autocorrelation function 



C vv {n) 



(2q - l) n dq = 



|(1 + (-!)") 
n + 1 



E„ 



n+1 



(9) 

where E m — 1 for m even and E m = for m odd. 

From the autocorrelation function, we can calculate 
the mean squared displacement by the general relation 

+n 

Ax2(n)=5x 2 ^2 C vv {m){n-\m\), (10) 



Ax 2 (n)/5x 2 = n + 2n^2 CU(m-) ~ 2 E mC vv (m). 

(11) 



m—1 

In our special case, we obtain 



rn — 1 



Ax 2 (n)/6x 2 =n + 2n V" 



E„ 



m + 1 



E„ 



m + 1 



The first sum can be rewritten as 



A n = y ^"v= y 

m + i 

m=l,2,3,... m=2,4,6,. 



1 

m + 1 



(12) 
(13) 



Assuming that the time step n of interest is even, we 
obtain 



E n 



A n = y 

^ m + 1 

m=l,2,3,... fc=l,2,3,... 



n/2 

E 



1 

2k + 1 



T - n + 3 . _ , 3 . 

*(— )-n- 2 ) 



(14) 



where is the digamma function. The second sum 

can be rewritten as follows: 



E 



m E„ 



E 



771+1 L — ' 777 +1 

777 = 1,2,3,... 771 = 2,4,6,... 



Tl/2 

E 

fc=l,2,3,. 



2A- 



n 
2 



71/2 

E 



1 



fc=l,2,3,. 



2k + 1 



2k + 1 
A n (15) 



Inserting the sums into Eq.(12), one obtains 
Ax 2 {n)/Sx 2 = 77 + 2nA n - 2 ( - - A n 



= (77 + \)2A n 
= (77+ 1) 



T , 77 + 3 . 3 . 



(16) 



The case for odd 77 can be computed in a similar way. 
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